Hall conductance from Berry curvature in carbon nanotubes 
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We analytically show that a gap is induced around the Dirac point in the electronic spectrum 
of a previously metallic nanotube, in the presence of electric and magnetic fields perpendicular to 
the tube axis. For realistic values of the fields, a gap of at least a few meV can appear. Despite 
the quasi-one dimensional nature of the system, the gapped state is associated with a non-zero 
topological invariant and supports a Hall effect. This is revealed when the flux through the tube is 
varied by one flux quantum, which leads to exactly one electron per spin being transported between 
the ends of the tube. 

PACS numbers: 73.22.-f; 73.43.-f 



The behavior of electrons in carbon nanotubes un- 
der the influence of magnetic fields has received both 
theoretical-"^^ and experimental^ attention. When a 
strong magnetic field is applied perpendicular to the tube 
axis, the Fermi velocity of states around the Dirac points 
is reduced. This "flattening out" of the dispersion re- 
lation is reminiscent of the dispersionless Landau levels 
in a 2D sample exposed to a constant magnetic fleld. It 
is therefore natural to ask whether a carbon nanotube 
in a perpendicular magnetic fleld supports a Hall effect. 
One type of Hall effect has already been pointed out 
theoreticalljj^. In a strong magnetic field, there are states 
localized on the tube flanks (i.e. the parts of the tube 
wall that are tangential to the magnetic field). States 
on opposite flanks carry current in opposite directions. 
When the Fermi energy lies among these states, an elec- 
tric field perpendicular to both the magnetic field and the 
tube axis induces a current along the axis of the tube. 

In this letter we study a different type of Hall effect. In 
order to observe it, a gapped state is induced by means of 
an electric field. The mechanism that is responsible for 
the gap is the same as in Ref.[l. When the Fermi energy 
is in this gap, a circumferential electric field induces a 
current along the tube axis. At the heart of the effect is 
the anomalous velocity of Karplus and Luttinger^ and the 
Berry curvature of the band structure of Dirac fermions^. 
To make a link with the two dimensional topological in- 
variant that appear in this context, a further magnetic 
field, parallel to the tube axis is introduced. The result- 
ing flux through the tube plays the role of crystal momen- 
tum in the transverse direction in this one-dimensional 
systerrt^. 

Our main results are the following. Magnetic and elec- 
tric flelds perpendicular to the tube axis produce a gap 
around the Dirac point in the electronic spectrum (cf. 
Eq.[7a|). For instance, a magnetic fleld of SOT together 
with an electric field of 5 V//zm produces a gap of 3.7 me V 
in a tube with radius 2.5 nm. We derive an expression for 
the Hall conductance when the Fermi energy is in the gap 
(Ea.[T2|l. Finally, we show that varying the fiux through 
the tube by one fiux quantum (j)Q — h/e, leads to exactly 
one electron per spin being transported between the ends 



of the tube. The argument is the reverse of Laughlin's 
argument for the quantization of the Hall effect 

These effects occur when electrons are confined to a 
single graphene sheet that is rolled up into a cylinder. 
We would like to vary the flux through the cylinder by 
at least one flux quantum. In a 20 T field parallel to 
the cylinder axis, this requires a radius of 8 nm or more. 
Single wall tubes found to date have radii < 3 nir>ii and 
are therefore too narrow to contain a quantum of flux. 
Multi-wall carbon nanotubes on the other hand have 
outer radii of several tens of nanometers. There is ex- 
perimental evidence^ d^'^'^ that at low temperatures, and 
small bias voltages, transport in a multi-wall tube is con- 
flned to the outermost wall. We therefore expect that 
effects considered in this letter can be observed in multi- 
wall nanotubes. 

However, since a multi-wall tube is not the cleanest 
system to display the effect, we also mention another 
possible experimental method to obtain the system we 
consider. The idea is to manufacture an ultra- wide single 
wall tube starting from a multi-wall tube and manually 
removing all the inner cylinders. The removal of inner 
cylinders from a ^ 1 /im section of multi-wall tube has 
already been achieved experimentally^"*. All experiments 
that we know of, report the removal of inner cylinders 
to leave behind a casing of several outer cylinders. None 
the less, we see no fundamental obstacles to producing 
ultra-wide single wall tubes by this method. 

The system we consider is depicted in Fig. ([T]). In order 
to have a tractable problem, we work in the long wave- 
length limit, in which it is assumed that states contribut- 
ing to the Hall conductance have energies small compared 
to the hopping amplitude t ~ 2.7 eV. For the present we 
consider a single spin species. The effect of including spin 
is simply to produce a Zeeman splitting in the single par- 
ticle spectrum. We will include this effect eventually. 

According to Ref.[lB, electrons in the tube are de- 
scribed by the Dirac equation 

[vcr-{-ihdr + eA) + U]'<i' = e'<i!, (1) 

subject to the boundary condition 'i>{r + 2ttRx) = 
gi(e-27i-iyr/3)^^^-) where R is the tube radius. We have 
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FIG. 1: The system consists of a single graphene sheet 
wrapped into a cy Under of length A and radius R with R <^ A. 
A magnetic field B and electric field E that are parallel to 
each other, are applied perpendicular to the cylinder axis. 
There is an additional magnetic field parallel to the cylinder 
axis, such that the magnetic flux through the cylinder cross- 
section is $. 



chosen x to point along the circumference and y to point 
along the tube axis. The index v is determined by the 
tube chirality and assumes values of or ±1. (See Ref.fisl 
for the precise definition.) When 1^ = the tube is metal- 
lic, whereas it is semi-conducting if = ±1. The in- 
dex r distinguishes between valleys, with t = -1-1 for 
valley K and r = — 1 for valley K'. The phase 9 in 
the boundary condition takes into account the effect of 
the component of the magnetic field parallel to the tube 
axis. It is given by 6* = 27r$/0Oi where <& is the flux 
through the tube cross section. The vector potential 
A{x) = BoRshi{x/ R)y corresponds to the component of 
the magnetic field perpendicular to the tube wall. The 
scalar potential is U = Uq cos{x/R). 
We define k = k^x + kyij with 



k ™ 



e - 27ri/r/3 



(2) 



where m is an integer that labels transverse modes. 
We make the ansatz that eigenstates are of the form 
'^krii'r) = e''" ''V''fcj)(a;)/VA where A is the length of the 
tube. In the physically relevant regime where A is far 
larger than all other relevant length scales, results are in- 
dependent of the boundary condition in the y direction. 
For convenience, we therefore assume periodic bound- 
ary conditions in the y direction, so that ky is quantized 
as ky = 27rn/A with n an integer. The index 77 = ±1 
distinguishes between positive energy (conduction band) 
and negative energy (valence band) solutions. The wave 
function "0^^ is normalized to unity, and obeys periodic 
boundary conditions in the circumferential direction. It 
satisfies H{k)4'kri = E^feij where H{k) = Ho + Hi{k) and 



Ho = vcr ■ {—ihxdx + eA), 
Hi (fc) hvcr -k + U. 



(3a) 
(3b) 



To find the eigenstates of H{k) close to the Dirac point 
approximately, the strategy is the same as in Ref.d The 
two degenerate zero energy eigenstates of Hq are found 
exactly. Hi{k) is then treated as a perturbation. For this 
to be an accurate approximation, hvkx, hvky and Uq have 



to be small compared to the level spacing of Hq. In the 
limit of a weak magnetic field, i.e. Bq <^ fi/eR^ , the level 
spacing of Hq is given by Hv/R. In this case there is at 
most one transverse mode m, namely the one that mini- 
mizes \kx\, for which the first order perturbative result is 
accurate. (All other modes have ?iw|fca;| > Tiv/R.) In the 
opposite limit of a strong magnetic field, the level spacing 
of Hq is of the order of the magnetic energy v^/ehBo and 
the perturbation theory is accurate for many transverse 
modes m. 

The two zero energy eigenstates of Hq are 
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(4) 



where F = (R/lm)^ cos{x/R), and Z,„ = \/^/eBQ is 
the magnetic length associated with Bq. The normal- 
ization constant is given by = RIq{2R^ /1'^)/2it. 
To first order in k and C/q, the eigenstates of H{k) 
are ij^krjix) = xtr^f+{x) + Xkr^f-{x) where the spinors 

/ 4- 

Xkr) 



ixt^^Xk^V satisfy 



(hvcr.k + ^.az)xk± = ±?Xfer,- 



The energy of the state is 



e = ±e ±^/y? + {fiikf. 



(5) 



(6) 



The mass /i and renormalized Fermi velocity are respec- 
tively given by 



Ai = f/o/i(2i?7CJ/^(2i?74), 
V = 2itRN^v = i;//o(2i?V;^), 



(7a) 
(7b) 



where In{z) is a modified Bessel function. In Appendix 
[K\ we test the accuracy and regime of validity of Eqs. ([6]) 
and ([7]) by comparing to numerics. 

In our analysis, we have neglected the Zeeman split- 
ting. When the Zeeman splitting is included, the spec- 
trum is gapped when /i > hqBq, where /iq = 5.788 x 
10~^eV/T is the Bohr magenton. The gap is /AE — 
2(/i — fioBo)- For small enough Bq, fi ~ cUoR^Bq/Ti 
is linear in Bq. Writing Uq = cEqR where Eq is the 
electric field strength, the condition for a gapped spec- 
trum becomes Eq > TifiQ/e^R^. For a tube with radius 
R = 2.5 nm, a gap appears when Eq > 2.3 V//im. For 
instance, including the Zeeman splitting, an electric field 
of 5V/^m will produce a gap of 3.7 meV at Bq = 30T 
in a tube with radius R — 2.5 nm. Because the same Eq 
leads to a larger Uq for larger R, the Zeeman splitting 
becomes less relevant and at the same time attainable 
gaps become larger, in tubes with larger radii. 

The next step is to calculate the Hall conductance 
when the Fermi energy is in the gap. We use the ap- 
proximate expressions for the eigenstates and energies we 
obtained above. The Kubo formula for the contribution 
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of valley r to the Hall conductance reads 

^2 /-t„,\2 /-i fA pL /.A 



2e' [fivY 
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pL/ pl\ p Li pl\ 

/ dxi / dyi / dx2 / dy2 

Jo Jo Jo Jo 



Im{[*fei-(ri)t(Tj,«'fc,+ (ri)] [*fc,+(r2) V,*fc,_ (ra)] } 



(£1-62)2 

Using the perturbative solutions for ^k±, one finds 



(8) 



dx rfy^'fci,,i(r)t(T«'fc2^2(r) 
Jo Jo 



= 6k,M,2nRN'xl,^axk,r,,. (9) 

Substituting the explicit form of Xk± , the Kubo formula 
becomes 



k 



(10) 



The sum ;^ X^fe replaced by the integral ^ dfc to 
obtain 



(11) 



Summing the contributions from the two valleys together 
and substituting for kx from Eq. ([2]) leads to 
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72 + (27rm + e' + 27rj//3)2J' 

(12) 



with 7 = 2t: Rfi/hv. To obtain this result we tacitly made 
two assumptions. (1) We assumed that first order per- 
turbation theory in k and Uo is sufficiently accurate for 
all modes that contribute to ayx- (2) We assumed that 
modes included in the sum in Eq. (|12p but for which first 
order perturbation theory is not accurate, give a negli- 
gible contribution. Do these assumption place a restric- 
tion on the magnetic field strengths for which Eq. is 
valid? Comparison with numerical calculations (see be- 
low) indicate that the result is most accurate in the limits 
of small or large R/lm, and less accurate for intermediate 
values. 

In both the small and large 7 limits, further simplifi- 
cations are possible. For 7 ^ 27r, all but the term with 
the smallest denominator in Eq. (|12l) make negligible con- 
tributions to the sum. For metallic tubes, ly = and 
the Hall conductance as a function of 9 consists of a se- 
quence Lorentzian peaks at 9 = 27rm, to = 0, ±1, ±2, . . . 
with heights 2e'^ /jh and widths 7. For semi-conducting 
tubes, v — ±1 and there are two sequences of Lorentzian 
peaks in the Hall conductance, a.t 9 = 2-Km ± 27r/3, 



TO = 0, ±1, ±2, . . .. The peak heights are j^h and 
the peak widths 7. 

In the opposite limit of 7 ^ 27r, the sum in Ec^. 
can be converted into an integral, and we obtain 



7»27r _ 
'yx 



(13) 



The 9 independence of be understood as follows. 

Eigenstates close to the Dirac point are localized to the 
scale Im and therefore become insensitive to boundary 
conditions when Im < R- Since 9 only appears in the 
boundary condition, the gap is independent of 9. 

Eq. (IT51) shows that the Hall conductance ayx be- 
comes quantized when R/lm is large enough. Related 
to this, the average of <Jyx{9) over 9, i.e. {o'yx)g = 

d9 <7yx{9) /2-K is always quantized. Using the result 
of Eq. ([I^ for ayx{9), we obtain 



'yx/g 



(14) 



Although we have only derived this result in the small 
Uo limit, it holds beyond this regime. Indeed, it can 
be shown^ very generally that {ayx)g equals an integer 
multiple of /h. This quantization is topologically pro- 
tected. 

Now we compare our analytical results to numeri- 
cal results obtained using the nearest neighbor tight 
binding Hamiltonian of a zig-zag nanotube with 4000 
hexagons along the length of the tube, giving a tube 
length A = 0.85 /im. The computation time of the nu- 
merical algorithm limits us to considering fairly narrow 
tubes {R ~ 2nm). We therefore have to consider mag- 
netic fields that are unrealistically large in order to ob- 
tain values of i?/Z,n ^ 1. In physical realizations R will 
be larger, allowing the large R/lm limit to be reached for 
experimentally achievable values of Bq . 

In Fig. we plot ayx versus 9 for 9 between and 
TT. (The curves are symmetric about tt = so that 
we do not have to consider the interval (— tt, 0) sepa- 
rately.) In (a) and (b) are shown the analytical results 
of Eq. ([T^ (solid lines) and the small 7 approximation 
of a single Lorentzian peak (dashed lines) compared to 
numerical results (circles). In both cases Uq = 0.01 1 
and Bq = 0.005 h/ea^ ~ 53 T. In (a) the radius was 
R = 51a/27r so that L = 51a and ly ~ 51mod3 = 0. 
R/lm — 0.57 is small and thus we see a Lorentzian peak 
at 6* = 0. In (b) R = hOa/2-K so that L = 50a and 
V = 50niod3 = —1- R/lrn = 0.56 is again small and we see 
Lorentzian peaks at = ±27r/3. In both (a) and (b) we 
see good agreement between analytical and numerical re- 
sults. In (c) the magnetic field strength Bq — 0.015 ea? /% 
is in the intermediate regime {R/lm = 0.97). Here the 
analytical result does not reproduce the numerical result 
as well as before. This is due to the level-spacing of the 
zero order Hamiltonian being of the same order as the 
perturbation for some modes that contribute to the Hall 
conductance, (d) has Bq — 0.0'b ea^/h which is in the 
strong magnetic field regime {R/lm = 1-8). Eigenstates 
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FIG. 2: The Hall conductance ayx as a function of 9, the mag- 
netic flux through the tube, in units oiTi/e. The solid curves 
are the analytical result of Eq. (|12|) . The dashed curves for 
(a) and (b) show the single Lorentzian peak approximation, 
valid for small 7. The circles were calculated numerically 
using the nearest neighbor tight binding Hamiltonian for a 
zig-zag nanotube with 4000 hexagons along the length of the 
tube. Thus the tube length is A = 0.85 ^m. (a) and (b) show 
results for weak magnetic fields, i.e. R/lm small. For (a) a 
tube of radius R — 51a/2n, i.e. with 51 hexagons around the 
circumference, was used. For (b), (c) and (d) a tube with 
50 hexagons around the circumference was used. (This cor- 
responds to radii of 1.97 nm and 1.93 nm respectively.) In 
(a) and (b) Uo = 0.01 1 and the perpendicular magnetic field 
strength is Bo = 0.005ea'^ /h ~ 53T corresponding to a mag- 
netic length Im =— 3.5 nm. In (c) and (d) results are shown 
for intermediate and large R/lm- In (c) and (d), Uo = 0.007f. 
In (c) Bo = 0.015 ea^/?i corresponding to a magnetic length 
Im = 2.0 nm. In (d) Bo = 0.05ea^/?i corresponding to a 
magnetic length Z,„ = 1.1 nm. 



gaps of at least several meV can be induced. We derived 
an expression for the Hall conductance fEa. [T2)) when the 
Fermi energy is in the gap. The expression agrees well 
with numerical results from the nearest neighbor tight 
binding Hamiltonian. The non-zero Hall conductance 
leads to quantized transport. When the flux through the 
tube is varied by one flux quantum, exactly one electron 
per spin is transported between the ends of the tube. 

Appendix A: Gap behavior 

Here we investigate the behavior of the gap as a func- 
tion of the system parameters. We also test the accuracy 
and range of applicability of the perturbative results by 
comparing to results obtained by diagonalizing the near- 
est neighbor tight binding Hamiltonian for a zig-zag nan- 
otube numerically. 




— 7T n —7T 7T 

8 e 



are localized on the scale of Im and therefore insensitive 
to boundary conditions. The Hall conductance is quan- 
tized to e^/h (per spin) independent of 0. Note that the 
numerical results in (a) through (d) all obey the expected 
quantization of the average {<Jyx)g = e^//i (Eq. [T4)) . 

In order to observe the non-zero Hall conductance and 
the quantization of its average over 9, consider the fol- 
lowing experiment. Suppose the flux <J> through the tube 
is varied. This produces an emf —d^/dt = {h/e)d9/dt 
around the circumference of the tube. Owing to the non- 
zero Hall conductance, a current / is produced along 
the axis of the tube, such that / = {h/ e)(jyx{0)d9 / dt. 
The total charge transported through any cross-section 
of the tube as the flux is changed by one flux quantum, is 
Q = J^" dtl = (h/e) J^'' d9ayx{9). From Eq. ([Ill) then 
follows that Q = e, i.e. one electron (per spin) is trans- 
ported through the tube as the flux is changed by one 
flux quantum. This argument is Laughlin's argument for 
the quantization of the Hall effectiP applied in reverse. 

In conclusion, in the presence of parallel magnetic and 
electric fields that are perpendicular to the tube axis, a 
carbon nanotube that was previously metallic develops a 
gap. For realistic values of the fields and the tube radius. 



FIG. 3: Behavior of the gap. All plots are for an infinitely 
long zig-zag nanotube with radius 51 a. The solid lines are 
obtained from Eqs. ([6]) and (|7a|l . The circles are obtained by 
numerically diagonalizing the nearest neighbor tight-binding 
Hamiltonian of the tube, a) The gap AiJ versus the field 
strength Bo of the component of the magnetic field perpen- 
dicular to the tube axis. The electrostatic potential strength 
Uo ~ 0.01 1 is small enough for perturbation theory in Uo to be 
accurate. There is no magnetic fiux through the tube, so that 
9 = 0. b) The gap AE versus the scalar potential strength Uo- 
The magnetic field strength is Bo = 0.005 fi/ea^. There is no 
magnetic flux through the tube, so that 9 = 0. c) and d) The 
gap AE versus 9 (the flux though the tube in units of h/e). 
In both cases Uo = 0.007 1. In c) the magnetic fleld strength 
is Bo = 0.005 fi/ea^ so that the radius R is 0.56 times than 
the magnetic length Im and the gap is still strongly dependent 
on 6. In d) the magnetic fleld strength is Bo = 0.05 ti/ea^ so 
that the radius is 1.8 times larger than the magnetic length 
Im. This is sufficiently large for the gap to be independent of 
9. 

Results are plotted in Fig. ([3]). In (a) the gap Ai? is 
plotted as a function of Bq. The solid line shows the 
perturbative analytical result while the circles show the 
numerical result. The value of the scalar potential was 
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taken as Uq = 0.01 1. For small Bq the small parameter 
that controls the accuracy of the perturbation expansion 
is UoR/hv, which here has a value of 0.09. We see that 
this is sufficiently small for perturbative results to be 
accurate. In (b) the gap AE is plotted as a function 
of the potential Uq, at Bq ~ 0.005 h/ea^ and 2ttR = 
51 a. The analytical result co-incides with the numerical 
result for Uq up to ^ 0.04 i, corresponding to a value of 
UoR/hv = 0.37. In (c) and (d) the gap AE is plotted as 
a function of 9. In (c) a value of Bq — 0.005 h/ea^ was 
used corresponding to a ratio R/lm = 0.56. Since the 
magnetic length is longer than the radius, eigenstates 
are sensitive to boundary conditions and the gap has a 
strong 9 dependence. We also see that when 9 ^ zLtt, 



the first order perturbative result is not accurate. (This 
is because here the perturbation hvcr ■ k becomes large 
compared to the level spacing of Hq at small Bq.) In 
(d) a magnetic field Bq = 0.05 ?i/ea^ was used, which 
corresponds to a ratio R/lm = 1.8. Eigenstates close to 
the Dirac point are localized to the scale Z,„ and therefore 
become insensitive to boundary conditions when Im < R- 
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